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Abstract. This work explores the spectra of quantum graphs where the Schrodinger 
operator on the edges is equipped with a potential. The scattering approach, which 
(— I 1 was originally introduced for the potential free case, is extended to this case and 

O j used to derive a secular function whose zeros coincide with the eigenvalue spectrum. 

Exact trace formulas for both smooth and ^-potentials are derived, and an asymptotic 
semiclassical trace formula (for smooth potentials) is presented and discussed. 



1. Preliminaries and a summary of the main results 

A metric graph is a finite combinatorial graph where each edge is endowed with the 
natural metric and with a positive real length. A quantum graph is a metric graph 
together with a Schrodinger operator acting on it. The operator acts as second derivative 
plus lower order parts on the individual edges. Boundary conditions are imposed at 
the vertices to ensure that the operator is self adjoint. Quantum graphs where first 
introduced in the 1930s by Pauling to study the movement of electrons in molecules. 
They are a popular model for various physical processes involving wave propagation. 
The survey articles |GS06| . |Kuc0 8] and |Pos09j provide an excellent overview. 

One of the key properties that make quantum graphs such a useful model is the 
existence of an exact trace formula. The first trace formula for quantum graphs was 
proved in [Rot84j . It is valid for the standard Laplacian with Kirchhoff- Neumann 
boundary conditions at the vertices. Various generalizations have been shown, valid for 
Schrodinger operators with a magnetic field and much more general boundary conditions, 
(see |KS99bl IKPS07} IBE09] or |BE08] for a survey). However, all of these rely on the 
fact, that the eigenvalue equation can be solved explicitly on each of the individual 
edges, the eigenfunctions are linear combinations of counter-propagating free waves. 

Here, we consider the case of a more general Schrodinger operator A + w where 
w stands for the set of edge potentials, as will be detailed below. In this case the 
solutions on the edges are no longer known explicitly. However, Sturm-Liouville theory 
still guarantees the existence of two independent solutions on each edge. In the high 
energy limit these solutions converge to the freely propagating waves of the no potential 
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setting (see for example |PT87j ). Quantum graphs with a potential are also considered 
in |GS06j . |CS07] and [ HKT12] . Their methods are complementary to ours and the 
focus is on different aspects of the problem. 

Before stating our main results we provide some necessary definitions and describe 
the model. 

1.1. The model 

Let Q be a graph defined by a vertex set V and an edge set S. Their respective 
cardinalities will be denoted by V = \V\ and E = \£\. The natural Euclidean measure 
is attributed to the edges, and the length of an edge is denoted by L e , for e G £. The 
number of edges and vertices as well as the length of each edge are assumed to be 
finite. Distance along an edge can be measured from either of the connected vertices, 
and it is convenient to use a directed edge notation to specify the direction of increasing 
coordinate. Thus, given an edge e G £, one defines the two corresponding directed edges 
denoted by d and d. The set of directed edges will be denoted by T>. The initial and 
terminal vertices of a directed edge d G V will be denoted by i(d) and r(d), respectively. 
Clearly i(d) = r(d). The coordinate on the directed edge will be denoted by Xd and it 
assumes the value at i{d) and L d at r(<f). Thus, Xd = L d — x d . It is convenient to 
associate the length to the directed bonds, in which case L d = L^ = L e . 

The Schrodinger operator on each directed edge is defined as the sum of the 
(magnetic) Laplacian and a real continuous potential w d (x): 



Here A d are the position independent "magnetic potentials" associated to each edge 
with Ad = —Aj. For simplicity of exposition, it will be assumed that A d = on all 
edges, the generalization is straight forward. Denote the operator A w — A w . Clearly, 
w d {x) = w d (L e — x). The edge operator acts on the space of twice differentiable complex 
valued functions in the interval [0, L e ]. It is natural to define the edge potentials to be 
continuous on all the vertices: u^(0) = tt»<z'(0) for all d, d' : i(d) = i{d'). 

The graph Schrodinger operator acts on 2E dimensional vectors {ipd{.x)) d&v . To 
get a consistent description, they are invariant under a switch of the preferred direction 
on the edges, ipd(x) = ?Pd(Ld — x). The bond wave functions are continuous at the 
vertices and are doubly differentiable in (0, Ld). The Schrodinger operator is self adjoint 
if appropriate boundary conditions are chosen at the vertices. For the sake of clarity of 
notation, only the Kirchhoff-Neumann boundary conditions will be used: 



The formalism can be easily extended to include other types of boundary conditions as 




with d G V . 



(1) 




(2) 
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It is instructive to recall the main steps taken to derive the trace formula and point 
out the differences and similarities between the no potential done in |KS97j and 

[KS99bj and the potential case here. 

On every directed edge d the eigenvalue equation 

- ^'(x) + w d {x)Mx) = k 2 M%) (3) 
has two independent solutions, denoted by ip d an d V'd ■ Normalize them so that 

V£(o) = i (i£y(p) = T ik (4) 

This corresponds to the two solutions e T in the no potential setting, where the sign 
convention is made so that ip^ 1S the incoming wave and ip^ is the outgoing wave. The 
general solution is then given by 

tpd(x) = a d ipj(x) + b d ip^(x) (5) 

Finding an eigenfunction on the entire graph is now equivalent to finding a list of 
coefficients a d and b d . These coefficients have to satisfy two types of conditions: 

First, they have to be invariant under a switch of orientation on each edge. The 
general solution on the edge e can be rewritten in terms of r/>t and ipj, the two solutions 
on the same edge with reversed orientation. 

ip d (L e -x)= ifj^x) = a^t(x) + b#fjj{x) (6) 

Here x traverses the edge in the opposite direction. 

The resulting linear relation is captured in the edge transition matrix t^ e \k;w e ) 
defined by: 



a d 



t {e) (k;w e )[ b * ) (7) 



The matrix t^ (k; w e ) is unitary, it depends on the potential on the edge and on k. It 
will be derived explicitly in section |2j If there is no potential the edge transition matrix 
describes the phase shift acquired by free propagation along the edge, and t^ e \k;0) is 
a diagonal matrix with equal diagonal entries e lkLe . If there is a potential, there will 
be both transmission and reflection on the edge, and the edge transition matrix is not 
necessarily diagonal. 

Second, the eigenfunctions have to satisfy the boundary conditions at all the 
vertices, this step is equivalent to the no potential case and is captured in the vertex 
scattering matrices o v . 



»d— ^ °~dd> a d> (8) 
d':r(d')=v 

where t(d) = v. In physical terms, a v describes the scattering of waves at a vertex v 
as dictated by the vertex boundary conditions. For the Kirchhoff-Neumann boundary 
conditions it reads, 

^' = -V + d^)' for *-(<0 = *(<0 = « ( 9 ) 
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where deg(v ) is the degree (valency) of the vertex v and S is the Kronecker-5. Note that 
the scattering from an incoming edge to all outgoing edges is the same except for the 
scattering back to the same edge. 

The vertex scattering matrices are the building blocks for construction of the 
2E x 2E matrix 

^dd> = K{d%L(d)°M' where v = i{d) . (10) 

In a similar way, the graph transition matrix T(k; w) is constructed from the edge 
transition matrices 

T dd '(k; w) = (5^ d , + 6 d4 >)t d e J,(k; w e ), where e = (t(d), r(d)) . (11) 
1.2. The main results 

Putting all the above conditions and definitions together, gives rise to the secular 
function ^ w : C — > C The first main result of this work states its explicit form: 

Theorem 1.1. 

Let Q be a metric graph with the Schrddinger operator A w . Then there exists a secular 
function of the form 

( w (k) = (det S w (k))~^ det(Id 2E - S w (k)), where S w (k) = ET(k; w) , (12) 

whose zeros correspond to the eigenvalues of the operator including multiplicities (with 
the possible exception of the eigenvalue zero). 

Note that the secular function is of the same form as in the no potential case, see 
for example [KS99b] . The vertex scattering matrix £ is identical. The key difference is 
in the transition matrix T{k] w). If there is no potential this is a diagonal matrix which 
takes care of the phase shift accumulated when the wave propagates from one end of an 
edge to the other. In the present setting it is not necessarily diagonal, and it describes 
the transmission and the reflection induced by the edge potential. 

The Cauchy argument principle can be applied to count the zeros of the secular 
function ( Tl~2"I) and to derive a trace formula. However, before this is done, the concept 
of periodic orbits on the graph should be re-examined. 

In general, a periodic orbit is a closed oriented walk on the graph. It is completely 
determined by the topology (connectivity) of the graph. The periodic orbits on a 
potential free quantum graph coincide with the topological periodic orbits. However, in 
the presence of an edge potential, a classical particle or a wave is not only transmitted 
through the edge but can also be reflected back. Thus, a periodic orbit can include 
back-scattering from the edge potential - a dynamically induced variant which is not 
accounted for by the topology. In order to systematically include the possibility of 
potential induced reflections, the following scheme is proposed. It makes use of an 
auxiliary graph denoted by Q* . It is created from the original graph Q by inserting 
an auxiliary vertex (of degree 2) on each edge. The resulting graph is bipartite. Now, 
the reflection and transmission induced by the 2x2 transition matrices t^ e \k; w e ) can 
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be formally considered as the vertex scattering matrix of the corresponding auxiliary 
vertex. The set of periodic orbits in the trace formula is then the topological set of 
periodic orbits on Q*, which will be denoted by PO*. This is a new feature due to the 
edge potentials that does not occur in the free setting. 

The trace formula is a formal equality between distributions. To give it a meaning 
as a functional equality, a suitable space of test functions is required. Throughout, test 
functions are denoted by (p, it is assumed they are analytic on the real line, admit a 
holomorphic extension to a neighbourhood of the real line and are rapidly decreasing. 

The second main result of the present work is given in the following theorem. 

Theorem 1.2. 

Let Q be a metric graph with Schrddinger operator A w . Denote its positive eigenvalues 
by {k^} n . Then there exists K ^ such that 




for all test functions ip. The phase of det T(k; w) is denoted by Q w (k). The A p {k) are 
complex amplitudes associated to each periodic orbit. They are build from elements of 
the matrix S w (k) along the periodic orbit, the explicit formula is written down in ffify . 

The threshold K depends on the edge potentials and it is a necessary ingredient in 
the theory. Its origin will be explained both in the formal proof of the theorem and in 
the example where the trace formula for ^-potentials is derived. The trace formula ({TBI 
is exact and the coefficients Q w (k) and A p (k) depend implicitly on the potential. For 
high energies the eigenfunctions converge to the free solutions. Thus to leading order in 
k the trace formula above converges to the no potential case. Indeed, in the potential 
free case, K = 0, Qo(k) = 2kC, where £ = J2 e ee and the leading order part of 
A p {k) depends only on the vertex scattering matrix E and the edge lengths. Using the 
WKB-approximation it is also possible to compute the first order correction term. The 
resulting semi-classical trace formula will be discussed in the last section of the present 
work. 

The rest of the paper is arranged as follows. In section [2] the transition matrix 
for a single edge t(k; w) will be explicitly derived, and its relevant properties will be 
discussed. Using this information, the proof of the two main theorems will follow in 
section [33 The construction will be illustrated in section [4] by considering a simple 
system with a 5-potential on the edges. Finally, the asymptotic semi-classical trace 
formula will be derived in section [5j 

2. The transition matrix for a single edge 

Consider a single edge e of the quantum graph. As explained in section 11.11 the general 
solution of the eigenvalue equation can be expressed in two different ways corresponding 
to the different orientations of the edge. 

Mtf (x) + b d ^ d (x) = ajipi (L — x) + bjipj (L — x) (14) 
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Recall the definition of the edge transition matrix. 
Definition 2.1. 

The edge transition matrix t^ e \k;w e ) is defined by the equation 

= t^(k-w e ) ( b £ \ (15) 

Using the symbol ' to denote differentiation with respect to x, equation ( I14p implies 
^(0) = ipj(L) ^(0) = ~^{L) . (16) 

A direct computation then yields an explicit expression for the transition matrix. 
Lemma 2.2. 

The transition matrix is given by 

-wty+,i>2) -r d imp{L) - (^-)'(o)^t(D 

^+(0)(^)'(L) - (^+)'(0)^(L) -W(^,^ d ) 

where W(i/j^, ip^) = ^ti^l)' ~ (V'jOVj the Wronskian. 

In order to derive the trace formula the matrix ^(k^We) needs to be unitary on 
the real line. This is achieved through the choice of a suitable normalization of ip d an d 
ip d . The matrix t^ e '{k\ w e ) is unitary on the real line whenever 



^ + (0)=^(0) (18) 

which is satisfied with the normalization ^(0) = 1 and (f/>^)'(0) = ^fik. Note that 
the potential is real, so the system has two real solutions, hence ip d (x) = i/)~[(x) for all 
points on the interval. It also implies that W(ip^,ip^) G iR and \ip^\ = \i/J d \. 
The eigenvalues of t^ e \k; w e ) are 

For the derivation of the trace formula it is necessary to have control of the behaviour 
of t^ e \k] w e ) in the vicinity of the real k axis. This is done in using the following lemma. 

Lemma 2.3. 

For any given potential w e E L 2 ([0, L]) there exist a K such that the matrix t^ e '(k+ie; w e ) 
is sub-unitary for k > K and e > 0. 

For high energies the matrix t^(k;w e ) converges to the no potential case. 

&\k;w e )={ n eikL 1+0(0 (20) 
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Proof. Take k real. From standard properties of Sturm-Liouville problems, |PT87j , the 
two solutions ij) ± satisfy 



\(iP ± )'(x)T'ike Tlkx \ < 2\\w\\ L 2e lM t 2VI 



k 



(21) 



and 



~dk 



(x) — — i ip (x) (ijj + (x)ip (x) — ip + (x)ip (x)) dx 



±v 



dk 



(x) = —i / ip (x) {x)il)~ (x) — i/j + (x)(vjj~y(x)) dx 

Jo 



Define an error function E(x,k) via 

^+(x) = e~ ikx + E(x, k) 
then the bounds above imply 



rj,-(x) = e tkx + E(x,k) 



(22) 
(23) 



\E{x,k)\ < 



C 
~k 

C 



\E'(x,k)\ < C 



< C 



(24) 



for some constant C depending on the potential but independent of k and x G [0, L] . 
Plugging ( !23~1) in the transition matrix yields 

1 



t (e \k-w e ) 



-ikL 1 

2ik 



\r(E'(L,k)-ikE(L,k)) 



±-(E'(L,k)+ikE(L,k)) 



2ik 



hr[E'{L,k)-ikE(L,k) 



This implies 



&\k; 



e ikL Q 

e ikL 



Oik 



(25) 



(26) 



as claimed. The eigenvalues of t^ e \k;w e ) are given by 
l±±\E'(L,k)+ikE(L,k)\ 



A*l,2(fc) 



-ikL L 

2ik 



{E'{L,k) -ikE{L,k)) 



ikL 1 e lkL (E'(L, k) - ikE(L, k)) ± \E'(L, k) + ikE(L, k) \ 



2ik 



-ikL 1 

2ik 



1 - k (E'(L,k)-tkE(L,k)) 



Let 



1 e lkL (E'(L, k) - ikE(L, k)) ± \E'(L, k) + ikE(L, k)\ 



(27) 
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then 



< 



C 
J 



< 



c 
T 



(29) 



from the bounds on E(L,k) and E ; (L,k). Here C can be computed from C, it also 
depends on the potential but not on k or x. 
Thus 

d 



de 



\Hi, 2 (k + is) 



-2L + 0(k- 1 ) 



e=0 



which is smaller then zero for large enough values of k. 



(30) 
□ 



Remark 2.4. 

This lemma still holds in the case of a 5-potential as will be shown in section HI While 
one could use this lemma to derive an explicit bound for K the resulting bound would 
be very crude. 

In the case of no potential or a positive 5-potential one can choose K — 0. For a 
negative (^-potential the transition matrix t^ e \k; w e ) is not subunitary above the real axis 
for small values of k and an explicit value for the optimal value of K will be computed. 



3. The secular equation and trace formula for graphs with edge potentials 



We will now derive the secular function ( w for a quantum graph with potential. Each 
edge will be treated as a scatterer similar to the vertices, using the transition matrix 
derived in the previous section. The matrix £ controls the scattering at the vertices of 
the graph, and it is not affected by the introduction of the potential. On all the vertices 
Kirchhoff-Neumann boundary conditions are imposed. 

Using the 2E x 2E matrices £ and T(k;w) (llOflip we recall 

S w (k) := Tj ■ T(k; w) . (31) 

Now the conditions on the coefficients and bd can be written as 

f h\ ( h\ 



S w (k) 



\b D J 



(32) 



\b D J 



In other words, the eigenvalues of A w on the graph correspond to the values of k 
where S w (k) has eigenvalue 1, including multiplicities. The eigenvector defines the 
eigenfunction in terms of the ip ± . This argument does not hold for k = as the 
functions ip^ 1 are not linearly independent with our normalization. 

Note that S w (k) is unitary on the real axis and subunitary above the real axis for 
large enough values of k by lemma 12.31 and the fact that £ is unitary and independent 
of k. 

To prove theorem 11.11 recall the definition f|T2|) of the secular function 

C w (k) := (det S w (k))-* det(Id 2E - S w (k)) . (33) 
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This function is zero whenever S w (k) has eigenvalue 1, the prefactor is always nonzero 
because S w (k) is unitary. It is inserted to make ( w (k) real on the real axis. 

To prove theorem 11.21 one observes first that the function ( w can be extended 
holomorphically in k into a neighbourhood of the real line because the end values ^(L) 
depend holomorphically on k, [PT87] . Moreover, by the Schwarz reflection principle, 
the function ( w (k) satisfies 



C« (*) = Uk). (34) 

The spectral density is derived using the Cauchy argument principle. Let {k 2 } n 
denote the sequence of positive eigenvalues. In order to expand the spectral density into 
a sum of periodic orbits, the matrix S w (k) needs to be subunitary above the real axis. 
By lemma 12.31 this is only true for k sufficiently large, therefore only eigenvalues larger 
than some constant K are counted. Negative eigenvalues, ie imaginary values of k are 
also ignored. Assume K 2 is not an eigenvalue. Let 

C £)K ■= {z e C | -e < Im{z) < e, Re{z) > K} (35) 

and let <p be a test function, then 

00 1 f d 

Y] viK) = — lim / ip(k) — ln(( w (k))dk 

1 f°° d 

= lim Im / (p(k + i£) — \n(( w (k + ie))dk (36) 

7T e->0 J K dk 

where we used the Schwartz reflection principle. The integral over the interval (—e,e) 
vanishes in the limit e — > because K 2 is not an eigenvalue. 
Plugging in the definition of ( w (k) gives 

1 f°° d 
^ <p(k n ) = -^~ Im / (p(k) — \a(detS w (k))dk 

I POO Q 

lim im / f(k + is)— In det(Id 2 E — S w (k + ie))dk (37) 

7T e->o J K ok 

In the first integral the limit e — > commutes with the integral because S w (k) is unitary 
and thus the function inside the integral does not have any poles on the real axis. 

Let Q w (k) = \ lndet T(k; w), then Q w (k) is real on the real axis because T(k] w) is 
unitary. Then 

^Hdet S w (k)) = z^^- (38) 

because E is fc-independent. 

If K is sufficiently large, by lemma l2T3| the matrix S w is subunitary above the real 
axis, so it admits the expansion 

oo 1 

\n(det(Id 2E - S w (k + ie))) = - V - tv(S^(k + ie)) (39) 

L — ' n 

n=l 
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i r 00 00 1 

+ -lim/m / <z?(A; + ie)— ^ - ti(S"(k + ie))dk (40) 
7r e^o J K ok ^— ' n 

n=l 

The tr(5^) terms can be interpreted as a sum over periodic orbits by using the identity. 



y tr(ST(fc;^))" _ ^ 1 / £ \ ^ 



n 

n=l 




Each edge in the original quantum graph is seen as a scatterer with scattering matrix 
given by the transition matrix t^ e \k\w e ) corresponding to the potential on that edge. 
Thus the tr(S w (k) n ) terms are counting the periodic orbits of the bipartite graph Q* 
built from the original graph Q by inserting an extra vertex on each edge. Let PO* 
denote the set of periodic orbits of Q*, p G PO*. Let n p be the topological length of 
the periodic orbit p, that is the number of edges it traverses. Denote by p the primitive 
periodic orbit that p is a repetition of. Then 

a 00 1 00 1 f) 

/m lE- tr ( S »W n ) = E- E n p Im-Y[T M ,(k) (42) 

n=l n=l p£PO*,n p =n d£p 

Here d' is the bond that follows d in p. The coefficients Tdd'{k) are of two different types. 
If the vertex between them is a vertex v of the original graph then Tdd'(k) = o~ d v d ] is the 
vertex scattering coefficient at that vertex from the S matrix, (which is independent of 
k). If the vertex between them corresponds to one of the edges e of the original graph, 
then Tdd'{k) = t d e J,(k) from the corresponding edge transition matrix. 

The factor np comes from the fact that each periodic orbit is counted once for each 
starting vertex. 

Theorem 11.21 follows once we identify 

A p (k) = ^Im-^T\Tdd>(k) (43) 
n v ok xx 

p dep 

The trace formula can also be written in distribution form 

A w4^ + ;E^) <«) 



^— ' 27T Ok 7T 

k n >K p£PO* 



4. Example - edges dressed with ^-potentials 

To illustrate the formal results, we discuss here a graph dressed with (^-potential on the 
edges. The solutions ip^ are explicit so one can write down an explicit trace formula. 
This example shows clearly how an edge with a potential acts as a scatterer. Note 
that a ^-potential can also be modelled by introducing a vertex with suitable boundary 
conditions, there is an exact trace formula for these cases, sec [KS99b] or |BE09j . The 
results here recover the ones in these two papers. 
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First consider a single edge parametrized as [0, L], look at 

- tf'(x) + D5 X0 {x)^{x) = k 2 ij(x) (45) 
This implies that ip' has a jump discontinuity at x$ namely 

lim^'^o + e) - i)\x Q - e)) = D^(x ) (46) 

Here D is some real parameter that measures the magnitude and x G [0, L] is the 
location of the (^-potential. 

Assume for now that k is real, so only positive eigenvalues are considered. The two 
individual solutions are then of the form 

~ ikx x < x 

ib + (x) — I n n 

+ x>Xtl 

e tkx x < x 

^-f x ) = J D D (47) 

;i + ^)e--^e— e"- x > x 

where the coefficients in the linear combination where computed from the continuity of 
the solutions and the jump in their derivatives. 

The Wronskian is then W(ip + , = 2ik on the entire edge. In terms of the edge 
d with reverse orientation the 5-potential is situated at L — x , so the transition matrix 
is 



\ ( \ _D_g— 2ifexogiA;L 

z,Ud X0 ) = —— 

with eigenvalues 



t(k,D6 X0 ) = {i _^ )e _ ik I n_ e2tkX0e _ lkL «* 1 I (18) 



»i(k) = e ikL KM = f^f^ (49) 

Lemma 4.1. 

If the 5-potential is positive, D ^ ; the matrix t{k + ie\ D5 X0 ) is subunitary for all real 
k. If D < it is only subunitary for sufficiently large values of k, namely for 

m 

Proof. Let e > then 

fi^k + ie) = e ikL e' eL (51) 

so 

\ l i 1 {h + te)\<\n 1 (k)\ = \ (52) 
For the second eigenvalue 
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thus 



which means 



^- \fx 2 {k + ie)\ 2 
as 



s=o D 2 + 4fc* 



S/} 2L (54) 



M*+*OI<i * l> JdWT¥ (55) 

This is trivially satisfied if D ^ but not true in general if D < 0. □ 
Remark 4.2. 

The matrix t(k; D8 X0 ) is only subunitary above the real axis for large enough values of 
k. This is the same situation as in lemma 1231 for potentials in L 2 ([0,L]). 

The trace formula contains terms of the form Jrlndett(fc; D5 Xo ) which can be 
explicitly computed using 



2ik + D j kL _ 2iarctan(#)+iA:L 

so that 



^k)=e ikL /x 2 (fc) = - ~ e ikL = e - 2iarctan (^)+' fcL (56) 



I|lndet t (^) = 2L + i ^ (57) 
Corollary 4.3. 

The trace formula with a 5-potential of strength D e on each edge reads 

co r E 

E W*) = 7 + ;EiBX75i + ; E Mk) (58) 



7T 7T ^ Ak 2 + D 2 7T 

k n >K e=l e pePO* 



coefficients 



Ap(k) = —Im^-Y[r ddl (k) (59) 



and the Tdd'{k) from the matrix ST(fc; D5 XQ ) with t^(k;D e 5 XQ ) computed above in 
equation\J^ Here K has to be chosen so that the inequality in lemma 4-1 is satisfied on 
all edges. 



5. The high-energy limit 

Consider the high-energy limit, that is, assume that k is real and k 2 » \\w\\l°°(g)- The 
exact solutions will be approximated by the WKB-method. This requires some estimates 
on the WKB approximation so the main steps in its derivation will be recalled. The 
goal is to find an approximate solution to 

— ift"(x) + w(x)ip(x) = k 2 ?p(x) (60) 

for large values of k. Set 

p(x) := a/ k 2 — w{x) s(x) := / p{x')dx' (61) 

Jo 
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and then use the substitution 

ip(x) =:p- 1/2 (x)r](s(x)) (62) 
This is well defined as p(x) > by the assumptions. The equation then transforms to 

d 2 (\ 5 \ 

— 7](s) + rj(s) = - ( -w"(x(s))p~ A (x(s)) + —w'(x(s)) 2 p~ 6 (x(s)) j rj(s) (63) 

Let x( x ) '■= \ w " '(x)p~ 4 (x) + -^w'(x) 2 p~ 6 (x) to simplify notation. Equation (163]) is 
solved iteratively by a power series r](s) := J2j =0 Vj( s )- The two initial solutions are 
770 (s) := e ±ts and the general term is the solution of 

d 2 

g^Vj(s) + r]j(s) = -x{x(s))Vj-i(s) (64) 

This is a second order linear differential equation which can be solved explicitly through 
variation of constants. 

Vj (s) = -e is [ S e~ 2is ' I" e is " X (x(s"))7 lj ^(s")ds"ds' (65) 
Jo Jo 

This process converges if x is sufficiently small, which is satisfied for k sufficiently 
large. The WKB-approximation is simply 770, written in the x variable this yields the 
familiar 



where the constant y / p(0) term was inserted for normalization. 
Lemma 5.1. 

On each edge the WKB-solutions differ from the exact solutions by an error term that 
decays like k~ 2 for k sufficiently large. 

C 

\tI>Hx) -^wkb( x )\ < p 

\m'{x)-{^ KB )'{x)\ <j (67) 
with C independent of k and x G [0, L]. 

Proof. For k large, p(x) ~ k so s(x) ~ k and x( x ) ~ k~ A . The iteration then gives 
r)j(x) ~ k~ 2 i and (r]j)'(x) ~ k~ 2 i +1 so the results follow from 

3=1 

□ 

On top of approximating the exact solution by the WKB-solution, the derivatives 
of the WKB-solution will be approximated as follows. 

^wkb)\ x ) = Tik^ KB (x) + Oik' 1 ) (69) 
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Thus 

V&jcb(0) = 1 (tPwkb)'(Q) = Tik + OCA;- 1 ) (70) 

Both layers of approximation have an error term that is two orders smaller than the 
main term. 

The transition matrix for the WKB-approximation can now be computed explicitly. 

t(k; w) 
= t WKB (k-w) + 0(k- 2 ) 





e * So p{x)dx 



e iJ L p(x)dx 




0(k~ 



+ 0(k- 2 ) (71) 



This information will now be used to study the behaviour of the trace formula in 
the high energy limit. Notice that the function s defined in (IBTI) . now seen as a function 
of k, corresponds to the classical action along an edge of the graph: 



Sd(k) = sj(k) := / \/k 2 - w d (x)dx (72) 
Jo 

Recall that Q w (k) = \ lndetT(fc; w), which can be approximated to 0{kr 2 ) by 

Q w (k) = Q% KB (k) + 0{k- 2 ) = s d (k) + 0{k- 2 ) (73) 

dev 

Note that the trace formula as stated in equation ( [TBI contains the derivative of 
the Q w (k), which in the WKB approximation reads 

i e -" w =gr^^ + i° ( ^- <74) 

Each of the integrals above can be identified as the time it takes a particle to traverse 
the edge (in the system of units used here, "time" is expressed in units of length). In 
the physics literature, the derivative of Q w is known as the Wigner delay time which 
will be denoted here by T(k). The result above can be interpreted as the Wigner delay 
time for a graph. In the high energy limit, it approaches the total length of the graph. 

An important consequence of the WKB approximation to the edge transition matrix 
(171]) is that reflections by the edge potential can be neglected to order p . This fact 
allows us to neglect the contribution to the trace formula from all the periodic orbits in 
the graph that involve a reflection on an edge. In other words, only the periodic orbits 
from the original graph Q survive. Their weights A p (k) split into a vertex scattering 
part and the contribution from the edges. 

A p (k)=< KB (k) + ^0(k- 2 ). (75) 

and 




Aj KB {k) = ?*Im TT a v M , TT e l /o d >/* 5 =^)*» . ( 76 ) 
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Here the a dd , = —5bb' + giffa) are fr° m the Kirchhoff-Neumann boundary conditions as 
before. To cast the trace formula in a form similar to the semi-classical trace formula 
|Gut90j . it is convenient to define the following quantities for each periodic orbit of 
topological length n p which consists of r p repetitions of a primitive periodic orbit of 
length rip: 

Let Ap = | n^ep a dd' I an d let vp be the number of back scatters (where a dd , < 0) . 
Let Sp(k) = J2dep s d(k) which is the classical action along the periodic orbit, and let 
Tp(k) = -^Sp(k) denote the classical period. Then, 



as the classical probability to survive along the period orbit, and vp the index replacing 
the Maslov index. Again to zeroth order this simplifies to the no potential situation. 
Finally, 



where E(k) is the error term, it decays as E(k) = 0(k 2 ). 
Remark 5.2. 

Note that the term involving the Euler characteristic from the classic trace formula is 
not visible because the formula here is only valid for large values of k and this term 
would appear as a 5 distribution located at k = 0. 
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